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ABSTRACT 

We present pre-perihelion infrared 8 to 31 /irn spectrophotometric and 
imaging observations of comet C/2012 K1 (Pan-STARRS), a dynamically new 
Oort Cloud comet, conducted with NASA’s Stratospheric Observatory for 
Infrared Astronomy (SOFIA) facility (+FORCAST) in 2014 June. As a “new” 
comet (first inner solar system passage), the coma grain population may be 
extremely pristine, unencumbered by a rime and insufficiently irradiated by 
the Sun to carbonize its surface organics. The comet exhibited a weak 10 /irn 
silicate feature ~ 1.18 ± 0.03 above the underlying best-fit 215.32 ± 0.95 K 
continuum blackbody. Thermal modeling of the observed spectral energy 
distribution indicates that the coma grains are fractally solid with a porosity 
factor D — 3 and the peak in the grain size distribution, a pea k = 0.6 /an, 
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large. The sub-micron coma grains are dominated by amorphous carbon, with 
a silicate-to-carbon ratio of 0.801 q2q. The silicate crystalline mass fraction is 
0.201 q^o, similar to with other dynamically new comets exhibiting weak 10 //rn 
silicate features. The bolometric dust albedo of the coma dust is 0.14 ± 0.01 at 
a phase angle of 34.76°, and the average dust production rate, corrected to zero 
phase, at the epoch of our observations was Afp ~ 5340 cm. 


Subject headings: comets: general - comets: individual (C/2012 K1 Pan- 
STARRS) - ISM: dust 


1. INTRODUCTION 


Solar System formation was an engine that simultaneously preserves and transforms 
interstellar medium (ISM) ices, organics, and dust grains into cometesimals, planetesimals 
and, ultimately, planets. Observing and modeling the properties of small, primitive bodies 
in the solar system whose origins lie beyond the water frost line (> 5 AU) provides critical 
insight into the formation of Solar System solids and establishes observatio n constra ints for 
planetary system formation invoking migration - th e ‘Grand Tack ’ epoc h flWalsh et al.ll201ll) . 


followed by the ‘Nice Model’ events flLevisonl 12009 : Gomes et al . 2005|). The characteristics 
of comet dust can provide evidence to validate the new, emerging picture of small body 
populations - including comet families - resulting from planetary migration in the early 
Solar System. 

Inside cometary nuclei, the bulk of the dust likely has been preserved since formation of 
the nucleus. Comet grains (and ices) also trace the pre-accretion history of comet materials 
extant in the outer disk. Comet dust composition can be studied via Stardust samples, se¬ 
lected collections of Interplanetary Dust Particles (IDPs), and in situ analysis in comet flyby 
and/or rendezvous missions. Dust species that are best explained as products o f aq ueous 


alteration (e.g., magnitite, cubanite, p ossibly pentlandite) are rare (IStodolna et al.l 12012 


Berger et ali l2 01ll : IZolenskv et al . 2008! ) and corresponding altered silicates (e.g., phyllosili- 
cates, smectite) are missing suggesting that aqueous alteration in cometary nuclei is limited, 
is no t we ll rep resent ed in the Stardust samples, or that these minerals have exogenous ori¬ 
gins (IBrownleel 2014 ). Thus, the bulk of comet grain properties including dust size, porosity, 
and composition relate to grain formation, radia l mixing, and particle agglomeration in the 
proto-solar disk (for an extensive review see Bro wn lee 20141. However, opportunities to study 
actual samples of cometary dust are rare, motivating the need for telescopic remote sensing 
observations of dust whenever apparitions are accessible from terrestrial observatories. 
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In this paper we report our pre-perihelion (TP = 2014 Aug 27.65 UT) infrared 8 to 31 
micron s pectro photometric observations of comet C/2012 K1 (Pan-STARRS), a dynamically 
new (see Oortill950l . fo r a de finition based on orbital elements) Oort Cloud comet - (l/a or9 ) = 
42.9 x 10 -6 AU -1 f Williams! 2015 ) - conducted with NASA’s Stratospheric Observatory for 
Infrared Astronomy (SOFIA) facility during a series of four flights over the period from 2014 
June 04 to 13 UT. Contemporaneous optical imaging observations are also presented. 


2. OBSERVATIONS 


2.1. Ground-based Optical Imaging 


Comet C/2012 K1 (Pan-STARRS) was observed on 2012 June 01.22 UT and again on 
June 04.24 UT with the 2.3-m Bok Telescope at the Kitt Peak National Observatory. The 
comet was at heliocentric distance (r^) of 1.74 AU and 1.71 AU, a geocentric distance (A) 
of 1.66 AU and 1.69 AU, a phase angle of 34.62° and 34.76°, for each date respectively. 


The images were obtained with the 90Prime camera flWilliams et al.l 2004J ), a prime 
focus imager built for the Bok Telescope. At the time of observation, the 90Prime camera 
utilized a thinned back-illuminated CCD detector with 4064 x 4064 pixels with a pixel size 
of 15.0 fxi n. At prime focus the camera pixel scale is 0.45" which yields a field of view of 30.5 
x 30.5 square-arcmin. The instrument was equipped with Cousins/Bessel system broadband 
V and R filters. Multiple exposures (23 images in R band and 9 images in V band of 30 
seconds each) were obtained of the nucleus and coma of the comet with the telescope tracking 
at the non-sidereal rate corresponding to the predicted motion of the comet provided by JPL 
Horizon^/] in an airmass range of 1.40 to 1.74. 

All images were corrected for overscan, bias and flat-fielding with standard IRAF0 rou¬ 
tines. The data was photometrically calibrated using eight field stars of various spectral 
types with known V and R magnitud es selected from the Naval Observatory Merged As¬ 
trometric Dataset (NOMAD) catalog flZacharias et ah 2004) on the same CCD amplifier as 
the comet. The standard deviation of the photometric V and R zero points derived from 
the average of the field stars is of order 1% and no color corrections for spectral type were 
applied. The average nightly seeing was ~ 2.2" in both bands. A single 30 sec exposure in 


1 http: / /ssd. j pi. nasa.gov/horizons. cgi 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso¬ 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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the R band obtained on 2014 June 04.24 UT is shown in Fig. Q] 


2.2. SOFIA 


Mid-infrared (rnid-IR) spectrophotometric observations of comet C/2012 K1 (Pan- 


(FORCAST; Herter et al. 

2012) mounted at the Nasmvth focus of the 2.5-rn telescope of 

the SOFIA Observatory ( 

Youner et al. 2012: Gelirz et al. 

2009). The data were acquired 

over a series of four Lights, originating from Palmdale, CA at altitudes of ~ 11.89 km in 


2014 June, that were conducted as part of our SOFIA Cycle 2 programs to observe comets 
(P.I. Woodward, AORJDs 01_001 and 02_0002). Details of all SOFIA observations and the 
orbital parameters of comet C/2012 K1 (Pan-STARRS) at those epochs are summarized in 
Table [T] 


FORCAST is a dual-channel rnid-IR imager and grism spectrometer operating from 5 
to 40 /irn. Light is fed to two 256 x 256 pixel blocked-impurity-band (BIB) arrays, each 
with a plate scale of 0.768" per pixel and a distortion corrected field of view of 3.2' x 3.4'. 
The Short Wavelength Camera (SWC) covers the spectral region from 5 to 25 /irn, while 
the Long Wavelength Camera (LWC) operates at wavelengths from 25 to 40 /tm. Imaging 
data can be acquired in either dual channel mode (with some loss of throughput due to the 
dichroic) or single channel mode. 

Imaging observations of C/2012 K1 (Pan-STARRS) in three Liters were conducted on 
the first flight series, prior to three flights dedicated to spectroscopy. Spectroscopic obser¬ 
vations of the comet used two grisms, one in the SWC (Gill) and one in the LWC (G227), 
and the instrument was configured using a long-slit (4.7" x 191") which yields a spectral 
resolution R = A/AA ~ 140-300. The comet was imaged in the SWC using the F197 Li¬ 
ter to position the target in the slit. Both imaging and spectroscopic data were obtained 
using a 2-point chop/nod in the Nod-Match-Chop (C2N) mode with 45" chop and 90" nod 
amplitudes at angles of 30°/210° in the equatorial reference frame. 


All FORCAST raw image data p r oducts were processed u sing the FORCAST_REDUX 
Data Pipeline, vl.O.lbeta fcf.. IClarke. Vacca. fe Shupindl2014r) . which employed the reduc¬ 
tion packages FORCAST_FSPEXTOOL, vl.1.0, and FORCAST_DRIP, vl.1.0. Processing 
of the raw spectroscopic data was performed using the same packages, with the exception of 
FORCAST_DRIP, which utilized vl.0.4. Details of the FORCAST_REDUX Data Pipeline 
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can be found in the Guest Investigator Handbook for FORCAST Data Products, Rev. 


2.3. SOFIA Imagery and Photometry 


Aperture photometry of the SOFIA image data of comet C/2012 K1 (Pan-STARRS) 
was performed on the Level 3 pipeli ne coadded ( *.COA) data products using the Aperture 
Photometry Tool (APT v2.4.7; Laher et al. 20120 . At all FORCAST Liter wavelengths, the 
comet exhibited extended emission beyond the PSF of point sources observed with FOR¬ 
CAST under optimal telescope jitter performance^ The photometry was therefore conducted 
using a circular aperture centroided on the photocenter of the comet nucleus. We used an 
aperture of radius 13 pixels, corresponding to 9.984", with a background aperture annulus 
of inner radius 30 pixels (23.58") and outer radius of 60 pixels (47.16"). This aperture, 
which is ~ 3x the nominal point-source FWHM, encompassed the majority of the emis- 
sio n o f the comet an d coma. Sky-annulus median subtraction (ATP Model B as described 


m 


L aher et al . 2012) was used in the computation of the source intensity. The systematic 


source intensity uncertainty was computed using a depth of coverage value equivalent to the 
number of coadded image frames. The dominant source of overall uncertainty in the im¬ 
age photometry were image gradients due to imperfect atmospheric background subtraction. 
The calibration factors (and associated uncertainties) applied to the resultant aperture sums 
were included in the Level 3 data distribution and were derived from the weighted average 
of 3 calibrator observations of /? And (2 each) and a Boo (1 each). The resultant SOFIA 
photometry is presented in Table [2j 

Due to turbulence, telescope jitter, and differing chop-nod patterns, i.e., the chopping 
difference between beams and the nodding of the entire telescope field-of-view (for a for a 
discussion and il l ustration of this st andard infrared observing technique with SOFIA - see 
Temi et al.ll2014t lYoung et al.l 120121) executed in flight, the multi-filter imagery data could 
not be used to generate color temperature maps due to the unstable PSF. 

During flights primarily devoted to obtaining grisrn data f£ )2.4p . images of the comet 
where obtained through a single Liter at 19.7 micron. Figure [2] shows the 19.7 /im surface 
brightness distribution of comet C/2012 K1 (Pan-STARRS) observed on 2014 June 13.17 UT. 
The nucleus is unresolved and azimuthally symmetric with a radial proLle FWHM of ~ 1.01" 


3 https://www.sofia.usra.edu/Science/DataProducts/ 
FORC AST_GI_Handbook_RevAl .pdf 

4 see http://www.sofia.usra.edu/Science/ 
ObserversHandbook/FORCAST.html §5.1.2 
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and the coma is extended and diffuse. Low surface brightness emission extends in a vector 
direction commensurate with that expected for a dust tail. 


2.4. SOFIA Spectra 


Three temporally distinct spectra of comet C/2012 K1 (Pan-STARRS) were obtained 
in both grism over a series of flight sequences spanning 6 days (Table_[H). Many c omets 


exhibit temporal variability in the infrared over periods of hours (e.g., IWooden et al. 2004) 


to days at relatively similar heliocentri c distances due to coma jets related to nucleus activity 
and/or nucleus rotation period (e.g., iKeller et al.l 120071 : iGehrz et al.l Il995l ) that produces 
observable changes in the observed SEDs. Inter-comparison of each SED over this period 
showed no substantial changes in overall continuum flux densities nor spectral features to 
within the uncertainty per spectral resolution element. In addition, the 19.7 /irn aperture 
photometry suggests also that the level of coma emission did not markedly change (cf., 
Table [2]). Apparently, comet C/2012 K1 (Pan-STARRS) was fairly quiescent in its infrared 
behavior during this epoch given our signal-to-noise ratio and aperture size (12,300 km 
radius). Thus, the three independent spectra were summed together in pipeline processing 
to produce an average spectral energy distribution (SED). A 3-point unweighted rectangular 
smoothing function was applied to this average SED to increase the point-to-point signal-to- 
noise ratio of the data product used in our thermal model spectral decomposition analysis. 
The calibrated data products do exhibit a few artifacts near the edges of the 17-27 /irn 
spectral order where a few data points deviate upwards (near 17 /mi) or downwards (near 
27 /irn) from the apparent spectral trend. The spectra of comet C/2012 K1 (Pan-STARRS) 
are presented in Fig. [3l 


3. RESULTS 


Taxonomically comet C/2012 K1 (Pan-STARRS) is a member of the dynamical comet 
family deno ted as nea rly isotro pic comets (NICs), also commonly referred to as Oort Cloud 
comets (cf.. IDones et al. 20041). The interior composition of the ecliptic comets (ECs) and 
the NICs likely are preserved during their residence in the Scattered Disk and the Oort 
Cloud, but their surfaces are subject to various processing effects. Modeling the coma dust 
properties provides insight into the origin and evolution of dynamic comet families. 
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3.1. Thermal Modeling of the Coma SED 


Thermal modeling of the observed thermal infrared SED of comets obtained using re¬ 
mote sensing techniques enables derivation of coma dust grain properties. In particular, 
SOFIA (+FORCAST) provides spectroscopic coverage with the Gill grism to the region 
9-12 /tm which contains features from amorphous and crystalline silicates (e.g., 11.2 /mi) 
and organic species (e.g., PAHs). The G227 grism spans 17.6-27.7 /rm, encompassing dis¬ 
crete resonances from crystalline silicates as well as spectra l signatu res from carbonates and 
phyllosilicates, putatively argued to be extant in comets flLisse et al.l 12006! ). The SED slope 
at long thermal ( > 15 /tm) wavelengths provides constraints on the abundance of the larger 
grain population in the coma. Observations in these spectral regimes are key to ascertaining 
the origins of silicates within the solar protoplanetary disk, and placing early solar disk evo¬ 
lution within the context of oth er circumstellar disk s observed today t hrough comparison to 
model and laboratory data (cf., iLindsav et al.ll2013t Koike et al.l 120101) . 


Modeling the mid-IR SED of C/2012 K1 (Pan-STARRS) yields estimates of the coma 
grain properties. We constrain the grain parameters by chi-squared fitting thermal emis¬ 
sion models to the observed spectrum. The grain parameters included in the modeling are 
size distributions (n(a) da), porosity, the crystalline mass fraction (i.e., the fraction of the 
coma silicate grains that are crystalline), and relative material abundances. The dust tem¬ 
perature is calculated assuming thermal equilibrium of the grains; wherein the composition 
(mineralogy), size, and heliocentric distance determine the temperature of the grains. 

Comet grains are dominated in co mposition by a handful of silicate-type materials 
flHanner fe Zolenskvll2010l: IWoodenl 120081) : Mg-Fe olivine- and pyroxene-types in amorphous 
(glassy) forms and their crystalline Mg-end-members f orster i te (Mg ^SiOj) and enstatite 
(MgSiOa). Cometary aggregates also conta i n organics (ISandford et al.l 120061 1 or amorphous- 


carbon-like materials (IMatrait et al.l l2008t iFormenkoval 119991) that may be the glue that 


holds the amorpho us and crystalline m aterials together I Flvnn et al.ll2013 : C'iesla fc Sandford 
2012 ). Our model (IHarker et al.ll2002l and references therein) uses five materials: amorphous 
olivine and amorphous pyroxene with broad 10, 18, and 20 /tm emission features, amorphous 
carbon with featureless emission, and crystalline olivine (Mg-rich) and orthopryoxene with 
narrow peaks. Broad and narrow resonances near 10 and 20 /tm are modeled by warm chon- 
dritic (50% Fe; 50% Mg) amorphous silicates (i.e., glasses) and strong 11.25, 19.5, and 24 /tm 
narrow features from cooler Mg-rich crystalline silicate materials. 

The amorphous carbon component in our dust model is repr esentative of several key 
dust species - e.g., elemental carbon dust flFomenkova et al. 19941 ). an organic component 
with C =C bonds, identif i ed by XANES spectra n ear 285 eV that can include amorphous 
carbon flFlvnn et al.ll2013l. 120031: IWirick et al.ll2009lh and possibly other carbonaceous grains; 





























































however, overall model results do not depend on this degeneracy. We do not specifi¬ 
cally include Fe-Ni sulfides (such as pyrrhotite or troilite) in our models nor carbonates 


or phyllosilicate-rich mate rials 
ysis of Stardust samples ( Naka m ura-Messen ger et al 


The latter materials have not been d etecte c 

Wooden _ 20081: 


2011 


in track anal- 


Zolenskv et al. 


2008) nor are they un equivocally evident in remote sensing data ( Bursentova et al. 20121: 
Woodward et al. 2007). Phyllosilicates, specifically smectites including montmorillonite, 


chlorite, and serpentine, have 18-23 p m resonances that worsen spectral fitting of comet 
C/1995 01 (Hale-Bopp) (Wooden et al. 1999 ). Hybrid IDPs may contain up to 10% smectite 
(INakamura-Messenger et al.l l201ll) . Smectitie is spectrally distinguishable from amorphous 


anhydrous olivine-ty pe and amorphous pyroxene materials (INakamura-Messenger et al.ll2011 


Wooden et al. 19991) . yet it is not required for spectral decomposition. 


While FeS-type grains are pres ent in IDPs (IBradlev fe Pail 120001) . meteoritics samples. 


and c omets grains, such as Wild 2 ( Heck et al. 20121: Zolenskv et al. 


20081; Velbel & Harvey 


20071 ). our SOFIA spectra (Fig. |3]) do not exhibit the b r oad 23 //m spec t ral features often 


associated with fi ne-grained FeS (IBursentova et al.l 120121 : iMin et al.l 120051: iHonv et al.l 12002 


Keller et a l. 20021) . Larger FeS particles would be spectrally indistinguishable from larger 
amorphous carbon particles at mid- to far-IR wavelengths, yet robust optical constants 
spanning visible through the far-IR are lacking for FeS due to measurement challenges of an 
inherently extremely absorbing material (L. Keller, private comm.). Thus, thermal modeling 
of FeS gains is unc erta in. Stardust sam ples appear to be richer in FeS and poorer in car¬ 
bonaceous matter ( Joswiak et al.l 1201 2). so there is no basis as yet to make an assumption 
about the relative abundance of FeS and amorphous carbon-like materials in comet comae. 
Hence we presume that the majority of absorbing materials in cometary dust re-radiating 
the observed infrared S EP is dominated by olivine, pyroxene, and carbonac eous (amorphous 
carbon-like) materials (IBrownleell20l4 IWoodenll2008l: IZolenskv et al.l 120081) . This presump¬ 
tion provides a foundation for comparing compositional similarities and diversities of comet 
dust composition derived from thermal models. 


The best-fit chi-square model results are summarized in Table [3j The model fit to the 
observed grism spectra with the corresponding spectral decomposition of grain components 
is presented in Fig. SI Mineralogically, the grains in the coma of C/2012 K1 (Pan-STARRS) 
are dominated by amorphous materials, especially carbon. Our models produce a Hanner 
(modified-power law) differential grain size distribution (HGSD)@ peaking with grains of 
radii a pea k = 0.6 /irn, indicating relative moderately larger grains are present, and the grain 
power-law slope N = 3.4. In a HGSD the small radii grains at the peak of the grain size 


5 This power law (in grain radii, a) is defined as n{a) da = (1 — a 0 / a) M (a 0 / a ) N ; where a 0 = 0.1 /nn and 
a P eak = a 0 (N + M)/N. 
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distribution dominate the surface area and the flux density. 


Grains in comets are likely fractal porous aggregates Schulz et ah 0201511 . The grain 
porosity (P versus the dust radius a), parameterized by D, is defined as P = ( a/0.1 /im)( D_3 ) 
with D = 3 for solid and D = 2.5 for highly porous grains ([Woodward et al.l 20111) . Grains 
in the coma of C/2012 K1 (Pan-STARRS) also are solid (the fractal porosity parameter 
D = 3.0). Solid grains are not unusual, 65% of the Stardust tracks are carrot-shaped from 
solid terminal particles (IHorz et al.l 120061 ). The sub-micron sized silicate-to-carbon ratio 
derived from our models is 0.801 q2o- The uncertainty in the parameters derived from the 
thermal models are at the 95% confidence level. 


3.2. The 10 /tm Silicate Emission Feature 


The 10 /im silicate feature in comet C/2012 K1 (Pan-STARRS) is qu ite weak com¬ 


pared to comets like C/1995 01 (Hale-Bopp) or 17P/Holmes (e.g.. IWooden et al.l Il999 


Watanabe et al.ll20091) . Following ( Sitko et al 


20041) . at 10.5 /iin we find the silicate emission 


(defined as \Fio/F^ Unuum ] ) is 1.18 ±0.03 above a blackbody curve fit to the observed grism 
spectra continua longwards of 12.5 /im. The best-fit blackbody is 7% = 215.32±0.95 K (using 
Gaussian weighted errors) and color excess, defined as T^(fit)/ (278K r^ 0 ' 5 ) is = 0.992T0.004. 
The normalized {F\/F\ T ) SED in the region near the silicate feature at 10 /iin is presented 
in Fig. [5j Typically data near 8 /tm are used to establish the blue-continua (AA7.7 — 8.4 /am). 
Our estimate of the local 10 /irn continua may yield slightly lower temperatures than an 
estimate that included 8 /inr photometry. 

The large value of a pea k inferred from the thermal modeling of the observed SED of comet 
C/2012 K1 (Pan-STARRS) in 2014 June is commensurate with the weak 10 /im silicate 
feature. Smaller grains (a pea k~ 0.3 /iin) produce higher contrast silicate features. Grains 
of greater porosity also produce higher contrast silicate features in the 10 /im band. Long 
period NICs have ‘typical’ HGSD slopes of 3.4 < IV <3.7 and silicate-to-amorphous carbon 
ratios 1. The grain size distribution slope of C/2012 K1 (Pan-STARRS), N = 3.4, is 
not atypical. The preponderance of larger sub-micron grains ( a pea k = 0.6 /im) in the coma 
of comet C/2012 K1 (Pan-STARRS) results in cooler radiating dust that contributes to 
the ‘continuum’ under the 10 /im silicate feature and to the far-infrared flux density (see 
Fig. U). The sub-micron mass fraction is dominated by amorphous carbon grains. Amorphous 
carbon has a featureless emission spectrum that extends through the 10 /tin region, so low 
amorphous silicate-to-car bon ratio also can weaken the silicate feature strength ([Wooden 
2008 ; Wooden et al. 2004 ). 
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3.3. The Silicate Crystalline Mass Fraction 


The mass fraction of silicate sub-micron grains tha t are crystalline in comet comae is a 


keystone for models of early pla net-forming processes flBockelee-Morvan et al.l 120021; ICiesla 


20071: Hugh es fe Armitage 20-10). This fraction is defined as 


silicates _ 


fcryst 


E 


TTlcryst,x 


j (rry) I yy) silicates \ 

= l y' cryst,x I 1 amorphous,x) 


(i) 


where m x is the mass of species x. Crystalline species in comet grains provide a record of 
the high temperature process that formed dust in the inner disk of the solar system and the 
large scale mixing that transported these hot nebular pr oducts to the cold comet forming 
zones. Crystals, their composition ( e.g.. IWooden 2008 ) and shape ( Lindsay et al. 2013 ) 
trace inner-solar disk conditions (e.g., iQgliore et al.l 201 lh and offer a view into the earliest 
planet-forming processes that occurred in our early Solar System. 

Crystals from the inner disk wer e transporte d out to the comet-forming regime and 
mixed with “amorphous” silicates (cf. Cielsa 2 0111) . The “amorphous” silic ates are thought 
to be outer disk materials that probab l y wer e inherited from th e ISM (IBrownleel 12014 


Watson et al. 2009 


Kemper et al. 2004 . 2005 : Li fe Drainel 2001 ) in the infall ph ase of 


the disk. They have non-stoichiometric compositions (GEMS-like, iMatsuno et al.l 12012 


Bradley fe Pail l2004 iBradlev et al.l Il999l. and references therein) that include the compo¬ 
sitional ranges of olivine [(Mg^ , Fe(i_ J/ )) 2 Si0 4 ], with y ~ 0.5 for amorphous olivine, and 
[(Mg x , Fe(i_ 2 ,)) Si0 3 ] with x ~ 0.5 for amorphous pyroxene-type materials. Crystals are 
identified by narrow IR emission features (e.g., 11.2, 19, 23.5, 27.5, 33 /irn) superposed 
on an underlying thermal continuum in remote sensing spectra. Crystalline silicates have 
been detected using remote sensi ng techniques in the dust comae of all comet classes in¬ 


cluding C/1995 01 (Hale-Bopp) (IWooden et al 


Impact coma of 9P/Tempel 1 (IHarker et a l. 


1999 


2005: 


Harker e 


Lisse et a. 


the fragmentation outburst of 17P/Holmes (IReach et al.l 120101) . and several other comets 


2006 


; al. 2002 L 2004a) the Dee p 


Ha rke r et al. 2007). 


(IKellev fe Woodenl 120091 ; Woodward et ah 20111). Amorphous silicates are also detected in 
these comets as well. Crystalline silicates are found in abundance in t he Stardust samp les of 
81P/Wild 2; however, the amorphous grains are d ifficult to identify (llshii et al.l 120081 ) and 
may be limited to the smallest dust grains (IBrownleel 12 01419 . 


Crystalline silicates are rare in the ISM; however, they account f or A 2.2 % of the to¬ 


tal silicate component in the directi on of the Galactic: Center (IKemper et al.l 12004 120051) 
and <5% along other lines-of-sight (Li & Drainc 120011 ). Solar System crystalline sili cates 
detected in comets must be formed in the early stages of our disk’s evolution (Brownlee 
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2014) . Crystalline silicates require T > 1000K to form through either gas phase condensa¬ 


tion or annealing of amorphous (glassy) silicate grain s flWoodenl 12008 


Wooden et ah 2005 


Davoisne et al.ll2006t lHenningil2003t iFabian et al.ll2000h implying that the crystalline silicates 
must have been proce ssed in the disk near the young Sun or in sho cks out to a maximum 
distance of 3 to 5 AU ( Harker fe Desch 2002 : Wehrstedt fe Gail 2008 ) . Post-formation, they 
were transported radially outward into the comet-formation zones flCharnoz fe Morbidelli 
20071 ) - a process t hat is apparently ubiquitous in observations of external protoplanetary 
disks lOlofssson et all (120100 . Glassy silicate spherules (GEMS) and crystals are seen in ag¬ 
gregates in cometary IDPs. Large ‘terminal particle’ crystals and sub-micron cryst a ls (crys¬ 


tallites) are components of aggregate grains captured in Stardust samples (IBrownlee et al. 
20121. hood: IZolenskv et ajhoosl l200fih . 


Thus to first order, the diversity of comet dust properties reflects the temporal and 
radial gradients in our Solar System’s early history and similarities and differences in dust 
characteristics, including f cry st, may provide observational tests of of planetary migration 
models within the early solar system during the epoch of planet formation that resulted in 
a variety of small body dynamical populations. We find a that the silicate crystalline mass 
fraction in comet C/2012 K1 (Pan-STARRS) is f cry st — 0 . 20 + qiq . This range is similar to 
that found for comet C/20007 N3 (T ulin), f £rus t = 0.48±0.06 (derived from the mass fractions 
presented in Table 3 of Woodward et al . 1201 ll. and Eqn. Q] of H3.3I) which also exhibited a 
weak 10 /irn silicate emission feature. 


3.4. EC and NIC Dust Characteristics 


As a result of giant planet migration, some comet nuclei were dynamically scattered into 
the Oort cloud to be exposed to the Galactic environment, whereas those bodies comprising 
the bulk of the ECs population have nuclei exposed and processed (at depths ranging from 
mm to few cm) by solar insolation, space weathering, and heliocentric activity variations 
(sublimation of CO, CO 2 ; crystallization of water and other ices) which affects materials 
lofted into the comae. Although the interior compositions of ECs and NICs likely are pre¬ 
served, their surfaces have differing processing histories. 


Typically, active comets (arising from a population of NICs dynamically derived from 
the Oort Cloud and moving on long-period orbits) exhibit high contrast 10 /irn silicate 
features. In contrast, short-perio d ECs (i .e., Jupiter-family comets) have, on average, lower 
10 /iin silicate features strengths (jSitko et al.l 120041 ). and are t hought to have lower activities 
(cf., the active area and active fraction measurements of A’Hearn et al. 119951 ). For decades, 
the low-activity of ECs has been attributed to the accumulation of a rime of insulating 
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larger grains that were launched on non-escape orbits (jJew ittl 120071 ). Thermal models that 
fit observed infrared spectra of comets reveal that high contrast silicate features arise from 
comae having a preponderance of sub-micron grains (IHarker et al.ll2002[ lHanner et al.lll994i) . 
Comae without these sub-micron grains have weaker silicate features. In individual comets, 
variations in the silicate feature strength have been s een on short time sca l es corresponding to 


the aperture-crossing tim es of jets or coma features flWooden et al.ll2004 IHarker et al.ll2005 


20071 ; iGicquel et al . 2012). These variations are best explained by changes in the differential 
grain size (n(a) da) or fluctuations in the silicate-to-carbon grain ratio. 

Differences between EC and NIC coma grain populations may a rise from the surface 


layers EC nuclei being “processed” or weathered (e.g., Li et ah 20151 ). Processing of ECs 


surfaces may result from their frequent perihelion passages that decreases surface volatiles 
and small grains and leads to the creation of rimes and dust mantles. Evidence suggesting 
such processing occurs over millennia may be found in the analysis of material excavated from 


comet 9P/Tempel 1 by Deep Impact : t 

re dust grains in the ejecta were smaller than those 

in the ambient coma ( 

darker et al. 

2007 

) and the immediate comet surface contained a layer 

of carbon rich grains 

Sueita et al. 

200 

5) and a dust mantle comprised of compact 20 pm- 

sized dust aggregates 

Kobavashi et al. 

2013j. However, this conjecture is not definitive as it 


unknown whether or not the impact location reflects the global surf a.ce dust properties of the 
nucleus. In ECs, the coma 10 pm silicate feature strengths are low flKellev fe Woodenll20091 ) 
and the dust pro duction rates are modest. However, w hen EC nuclei have either fragmented 
(i.e., 73P/SW3, IHarker et al.l l201ll: ISitko et al.l 1201 if). e xplosively released materials from 
subsurface cavities (i.e., 17P/Holmes. iReach et al.l 120101) . or have had s ubsurface materials 
excavated from depth (i.e., the 9P/Tempel 1 Deep Impact encounter, Harker et ah 20051) 
the IR SEDs exhibit 10 pm relatively strong silicate feature emission ( > 1.2) arising from 
a population of sub-micron size silicate grain species. Whether or not the strong 10 pm 
silicate features arise from the release of sub-micron sized grains or the disruption of loose 
aggregates of fine particles (e.g., through gas-pressure disruption or impact fragmentation) 
is not known. Indeed the silicate feature in 9P/Tempel 1 changed from a EC-like spectra to 
NIC-like spectra immediately after Deep Impact event, returning to an EC-like state several 
tens of hours later (cf., Harker et al. 120051) . 


NICs are cano nically considered to be more pristine with higher surface volatile abun¬ 
dance (cf.. I Wooden 20081 ) - the effects of dwell time in the Galactic environment being more 
benign. Also, NICs are often considered a homologous population lacking significant nucleus 
evolution. Inner solar system apparitions of these comets frequently result in brilliant comae, 
with large dust production rates and pronounced silicate feature emission at IR wavelengths. 
It is not entirely clear whether the highly active nuclear regions of NICs can spawn small 
sub-micron grains responsible for the silicate feature emission, either by heritage or by frag- 
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mentation induced within the gas acceleration zone. However, whether comet evolution, 
such as processing in the Galactic environment, can be ignored when comparing the Oort 
cloud comet dust composition (including that expressed in f crys t ) is an open question. 


Table [4] present estimates of f cry st and select characteristics of the dust derived from 
thermal modeling of the mid-IR SEDs for a set of well-studied Oort cloud and “disrupted” 
Jupiter-family comets. The crystalline silicate fraction ranges appreciably, from ~ 10% to 
~ 80%. The compositional similarity suggests that Oort cloud and Jupiter-family comets 
have common origin sites within the earl y solar syste m (an argument that parallels that 
derived from volatile composition studies, A’ Hearn e t al. 2012), but the range of f cr yst values 
in Oort cloud comets suggest this class may be sampling a particular region that is not 
represented in the Jupiter-family members. This inference is intriguing; however, limited in 
robustness as any tentative conclusions are based on a limited sample size. Large sample 
sizes are required to substantiate or vitiate these trends. 


Comet C/2012 K1 (Pan-STARRS) and C/2007 N3 (Lulin) have modest mean values for 
fcryst (<48%), low silicate-to-carbon ratios, and grain size distributions that peak at large 
radii, >0.6 /irn. The 10 pm silicate feature is weak and/or absent in these NICs. Perhaps 
these bodies represent a population of more carbon dominated bodies, similar to the dark 
organic KBOs, whose surfaces are devoid of small grains. Indeed the low albedo (see f 33.6[) of 
C/2012 K1 (Pan-STARRS) and the dominance of amorphous carbon grain materials maybe 
providing clues. 


3.5. Dust Production Rates 


The radial profile of comet C/2012 K1 (Pan-STARRS) was plotted to assess the quality 
of the data for calculating a dust production rate near the epoch of our SOFIA observations. 
The radial profile of C/2012 K1 (Pan-STARRS) in the V band shows a deviation from the 
1/p profile flGehrz fe Nevl Il992l h suggesting contamination from gas su ch asCh (A = 0) 
band(s) near 5141 A. Strong C *2 emission is present in spectra flMcKav et al. 12 014 and also 
A. McKay, priv. comm.) contemporaneous with our optical imagery. We therefore only 
calculate the dust production in R band. The R band radial profile of C/2012 K1 (Pan- 
STARRS) is shown in Fig. [6] 

To estimate the rate of dust produc tion in c omet C/2012 K1 (Pan-STARRS), we utilize 
the Afp quantity introduced by A’Hearn et al. (11984T ). This quantity serves as a proxy for 
dust production and when the cometary coma is in steady state, the value for A(Q)fp is an 
aperture independent parameter, 
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A /p 2 y\ 2 1 rj 0.4(?77.come£ ^o) 

A(Q)fp = --- 

P 


cm) 


( 2 ) 


where A(0) is four times the geometric albedo at a phase angle 0, / is the filling factor of 
the coma, m comet is the measured cometary magnitude, m Q is the apparent solar magnitude, 
p is the linear radius of the aperture at the comet’s position (cm) and and A are the 
heliocentric and geocentric distances measured in AU and cm, respectively. To cor rect ou r 
comet me asurements for phase angle effects we applied the Halley-Marcus (HM) flMarcus 
2007a,b;; Schleicher et al. 1998) phase angle correction!®] We adopt an interpolated value of 
0.3864 (appropriate for the 2014 June 04.24 UT dataset) to normalize A(Q)fp to 0° phase 
angle. Table [5] reports values of Afp = [(A(0)/p/HM] at a selection of distances from the 
comet photocenter in the A-band. 

In addition to Afp , we also compute the efp parameter of comet C/2012 K1 (Pan- 
STARRS) based on our FORCAST broadband photometry (Table [ 2 ]). The efp parameter 


(defined by iKe llev et ah 20131 . Appendix A) can be considered to be the thermal emission 


corollary to the scattered-light-based Afp: 


A 2 F 

efp cm), (3) 

TTP By 

where e is the effective dust emissivity, F u is the flux density (Jy) of the comet within 
the aperture p, B u is the Planck function (Jy/sr) evaluated at the temperature T = 
Fcaie (278 K) r^ 0 ' 5 , where the scaling factor T sca i e = 0.99 based on the 215 K measured 
continuum temperature discussed in 1 13.21 Derived values of efp for comet C/2012 K1 (Pan- 
STARRS) are presented in Table [2l 


3.6. Coma Averaged Dust Albedo 

Dust albedo is a basic parameter characterizing the size distribution and physical prop- 
that is, surprisingly, infrequently measured. Following the convention of 
the bolometric albedo, (. A bo i ometric = (Energy scattered /Energy incident ) is 


erties of comet dust 


Gehrz fc Nev (119921 ) 


A 


bolometric — 


/(e) 

i + /(©)’ 


( 4 ) 


6 see http://asteroid.lowell.edu/comet/dustphase.html 
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where for comet dust the incident energy is the sum of the energy scattered by the coma 
plus the total energy of the coma’s thermal emission at an observed phase angle 0 (Sun- 
comet-observer angle). The term /(©) can be determined from fitting the observed spectral 
energy distribution of the comet with appropriate Planck blackbody functions in the infrared 
(thermal dust emission) and reflected solar spectra at optical (scattering) wavelengths 


/(©) = 


[AF 


Aj max, scattering 


[AF 


X\max,IR 


(5) 


where the [AFAjmaa; is the peak of the SEDs in the respective wavelength ranges. Lab ex¬ 
periments and theoretical calculations of the scattered light from particles indicate that the 
total brightness, color, polarization, and polarization color depend on the optical constants, 
particle size distribu t ion, structure, and porosity of the dust as well as the solar phase angle 
flLindsav et al.l 120131: lHadamcik et alJl2007t iKolokolova et al.l 120041) . The spectral shape of 
the IR thermal emission provides a direct link with the mineralogy and grain size. Both of 
these processes provide information on the size and composition of the dust. The scattered 
light and thermal emission are also connected to one another through the grain albedo, the 
ratio of the scattered light to the total incident radiation. Because light is not isotropically 
scattered by comet dust the measured albedo will depend not only on the composition and 
structure of the dust grains, but also on the phase angle (Sun-comet-observer angle) of the 
observations. 

The coma SED of comet C/2012 K1 (Pan-STARRS) was measured on 2014 June 04 UT 
using filter photometry at both mid-IR as well as the optical (scattered sunlight ) wav elengths. 
These data enable computation of the coma averaged bolometric albedo ( Gehrz fc Nevil l 992l ). 
Using the integrated flux densities in a circular aperture of radius 9.989", [AF A ] ma a;,/j? was 
derived from the SOFIA photometry by x-square fitting a blackbody to the mid-IR data 
using Gaussian weighted errors, resulting in T&& = 214.04 ± 14.94 K with a peak flux of 
2 . 02 / 0/5 x 10 ~ 16 W cm -2 . The [V] band photometry is contaminated by gas emission (]l )3.5p . 
However, the C/ bands fall outside the bandpass of the [R] filter and a 5800 K blackbody 
(the Sun) emission peaks near the [R] filter central wavelength (A c = 0.64 /mi) in AFa 
(W cm -2 ) space. Hence, [AFa] max,scattering = 3.33 ± 0.03 x 10 -17 W cm -2 derived the [R] 
band photometry measured in a circular aperture of radius 9.989" (Table [5]). The dust 
bolometric albedo of comet C/2012 K1 (Pan-STARRS) is 0.14 ± 0.01 at phase angle of 
34.76° (from Eqns. [4] and [5]). 


Kolokolova et al. (120041) reviewed published visual albedos of comets and found only 


eight 8 c omets have measur ed albedos (excluding comets Kohoutek an d Crommelin dis- 
cussed in lGehrz fe Nevlll992l) . and all were from the NIC dynamical class. iKellev fe Wooden 
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(120091 ) found only one EC with visual albedo, 21P/Giacobini-Zinner flPitti chova et al. 2008! ). 
Recently, the albedos of 73P/Schwas smann-Wachmann 3, 103P/Hartley 2, and C/2009 PI 
(Garradd) also have been measured (IMeech et all 1201 It ISitko et al.ll201ll. 120131) . Fig. [7] is a 
compilation of the the bolometric albedo data that exists on comets, including our determi¬ 
nation for comet C/2012 K1 (Pan-STARRS). There is considerable scatter for multi-epoch 
observations of individual comets. Such scatter arises from variations in activity of a comet 
at different epochs of observation. For example, comet C/1995 01 (Hale-Bopp), whose data 
are most scatt ered, had numerou s and fast changing mo rphological structures (jets, shells, 
envelopes, e.g., lHarker et al.lll997l; IWoodward et al.lll9981). All of these features were charac¬ 


terize d by differing size and particle composition (e.g.. iRodriguez et alJll997tlSchleicher et ah 


1997 ). Thus, the difference in the dust albedo for the same comet indicates variations in 
comet activity, specifically development of jets and other morphological features. 


The ensemble albedos compiled by iKolokolova. et al.l (120041 ) also shows a broad distri¬ 
bution of values for each phase angle. The causes for these latter albedo ranges and the 
scatter in multi-epoch observations of comets are unclear, but must reside in the physical 
properties of the comet particles, including the grain size distribution, porosit y, grai n struc¬ 


ture (i.e., prolate spheroids, crystalline needles, etc.), and composition (e.g.. iLindsav et al. 


20131 ). However, observations have not yet demonstrated to what extent grain structure or 
grain compositions are important. To assess these latter aspects, thermal emission models 
and albedo observations of a additional comets are needed. 


4. SUMMARY 

We discuss the pre-perihelion mid-infrared spectrophotometry and narrow band filter 
imagery obtained in 2014 June with FORCAST on the NASA SOFIA airborne platform 
of the dynamically new comet C/2012 K1 (Pan-STARRS) at a heliocentric distance of ~ 
1.70 AU. The spectral energy distribution of the comet at this epoch exhibits a 10 //in silicate 
feature, [Fio/F^ tinuum ] = 1.18 ± 0.03 above a blackbody curve ( T bb = 215.32 ± 0.95 K) fit 
to the spectra continua longwards of 12.5 /im which is quite weak compared to comets such 
as C/1999 01 (Hale-Bopp) or 17P/Holmes. The coma dust bolometric albedo, 0.14 ± 0.01, 
derived using contemporaneous optical imagery is similar to other comets at the observed 
phase angle (~ 35°), while the dust production rate ( Afp ) from scattered light observations 
is ~ 5340 cm. 

From the observed infrared spectral energy distribution, thermal modeling analysis is 
used to determine the physical characteristics of the coma dust population and deduce the 
silicate crystalline mass fraction (0.20lg;fo) and silicate-to-carbon dust ratio (0.80/qJo)- We 
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find that grains in the coma of C/2012 K1 (Pan-STARRS) are dominated by amorphous 
materials, especially carbon, and the differential grain size distribution peaks at radii of 
0.6 //in, the slope of the distribution N = 3.4, and the grains are solid, having a fractal 
porosity parameter D = 3.0. The bulk grain properties of comet C/2012 K1 (Pan-STARRS) 
are comparable to other Nearly Isotropic comets (NICs) with weak 10 /irn silicate features and 
similar in respect to coma grains seen in the small-set of Ecliptic family comets (ECs) that 
have fragmented, explosively released subsurface materials, or have had materials excavated 
from depth. 

SOFIA observations of comet C/2012 K1 (Pan-STARRS) and other future comets en¬ 
ables characterization grain properties in the NIC and EC dynamical families. These proper¬ 
ties, including dust size, porosity, and composition, relate to grain formation, radial mixing, 
and particle agglomeration in the proto-solar disk and provide insight to the evolution of the 
early solar system. As the number of well-studied comets increases at infrared wavelengths 
(from whence dust properties can be characterized), the fundamental differences between 
comets originating from different regions and times in the solar system may be eventually 
discerned. 
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Fig. 1.— The Bok i?,-band optical images of comet C/2012 K1 (Pan-STARRS) obtained 
on 2014 June 04.24 UT, with logarithmic greyscale color map. The vector indicating the 
direction of the comet’s motion and the vector indicating the direction toward the Sun are 
also provided. 
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Fig. 2.— The SOFIA 19.7 /im image of comet C/2012 K1 (Pan-STARRS) obtained on 
2014 June 13.17 UT, with logarithmic color map. The vector indicating the direction of the 
comet’s motion and the vector indicating the direction toward the Sun are also provided. 
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Fig. 3.— The SOFIA grism observations of comet C/2012 K1 (Pan-STARRS). Grism ob¬ 
servations from three flight series where combined to create a single average SED, and a 
three-point unweighted rectangular smoothing function was applied to the composite spec¬ 
tra to improve the signal-to-noise. The filled circles are the broadband photometric obser¬ 
vations obtained on 2015 June 04 UT when each of the Elters were observed on a single 
flight. The photometry is scaled by a factor of 0.75 (which is of the order of absolute photo¬ 
metric calibration uncertainty of FORCAST) to match the flux density of the grism spectra 
and demonstrate that the shape of the photometric SEDs mimics that of the spectroscopic 
spectral segments. The vertical grey bar indicates the region of terrestrial ozone absorption. 
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Fig. 4.— The spectral decomposition of the SOFIA grism observations (see Fig. [3]) of comet 
C/2012 K1 (Pan-STARRS) derived from thermal modeling. The filled circles are the broad¬ 
band photometric observations obtained on 2015 June 04 UT when each of the filters were 
observed on a single flight. The filter photometry is scaled by a factor of 0.75 (which is of 
the order of absolute photometric calibration uncertainty of FORCAST) to match the flux 
density of the grism spectra and demonstrate that the shape of the photometric SEDs mim¬ 
ics that of the spectroscopic spectral segments. The vertical grey bar indicates the region of 
terrestrial ozone absorption. 
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Fig. 5.— The observed SOFIA grism flux density of comet C/2012 K1 (Pan-STARRS) near 
the 10 /mi silicate emission feature divided by a ~ 215 K blackbody continuum (Fx/F\,t) to 
highlight the details of the 10 /irn silicate feature. The vertical grey bar indicates the region 
of terrestrial ozone absorption. 
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Fig. 6.— Azimuthally averaged radial profile fluxes as a function of linear radius as measured 
in the R band from the optical photocenter of comet C/2012 K1 (Pan-STARRS) obtained 
on 2014 June 04.24 UT. The solid blue line denotes a 1/p profile. 



30 



Phase Angle (degrees) 


Fig. 7.— The bolometric albedo as a function of phase angle for a sample of Nearly Isotropic 
Comets (NICs; filled orange circles) and Ecliptic Comets (ECs; filled p urple squares) derived 
from the literature following the prescription of flGehrz fe Nevlll992l ). Our measurement of 
the NIC C/2012 K1 (Pan-STARRS), 0.14±0.01 at a phase angle of 34.76°, is indicated by the 
red star. The phase angles for each comet are obtained from the JPL Horizons ephemerides. 






Table 1. SOFIA Observational Summary - Comet C/2012 K1 (PAN-STARRS) 


Grisrn Total 


Observation 

or 


On Src 




Tail 3, 

Tail 3 

Date 

Fltr 

Exp 

Integ 



Phase 

Gas 

Dust 

2014 UT 

InstCfg A c 

Time 

Time 

Th 

A 

Ang 

PS Ang 

PsAMV 

(dd-mm hr:min:s) 

(>m) 

(sec) 

(sec) 

(AU) 

(AU) 

(°) 

(°) 

(°) 


FOF176 

06-04T03:35:31 

Imaging Dual 

19.71 

29.5 

616.0 

1.708 

1.688 

34.76 

105.75 

78.33 

06-04T03:35:31 

Imaging Dual 

31.46 

29.5 

616.0 






06-04T04:16:38 

Imaging Dual 

11.09 

30.8 

216.0 






06-04T04:16:38 

Imaging Dual 

31.46 

30.8 

252.0 






FOF177 

06-06T04:15:14 

Imaging SWC 

19.71 

45.0 

198.0 

1.684 

1.711 

34.76 

104.16 

77.47 

06-06T04:29:07 

Grism LWC 

G227 

22.5 

1800.0 






06-06T04:51:42 

Grism SWC 

Gill 

24.0 

1920.0 






FOF178 

06-llT03:53:40 

Imaging SWC 

19.71 

42.2 

204.0 

1.628 

1.769 

34.46 

100.71 

75.73 

06-llT03:57:33 

Grism LWC 

G227 

23.0 

1748.0 






06-llT04:32:22 

Grism SWC 

Gill 

24.0 

3336.0 






FOF179 

06-13T04:04:51 

Imaging SWC 

19.71 

45.1 

330.0 

1.605 

1.793 

34.22 

99.62 

75.17 

06-13T04:23:04 

Grism LWC 

G227 

23.0 

1564.0 






06-13T04:46:21 

Grism SWC 

Gill 

24.0 

2016.0 







*Notes. Orbital elements derived from JPL Horizons, ssd.jpl.nasa.gov/horizons.cgi. 
a Vector direction measured CCW (eastward) from celestial north on the plane of the sky. 
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Table 2. SOFIA Aperture Photometry and e/p of Comet C/2012 K1 (PAN-STARRS) 


Observation 

Date 

UT 2014 

(dd-mm hr:min:s) 

InstCfg 

(Imaging) 

Fltr 

A c 

(pm) 

Flux 

Density 11 

(Jys) 

XF x 

(xl0“ 16 W cm" 2 ) 

e/d 

(cm) 

FOF176 






06-04T03:35:31 

DUAL 

19.71 

15.102 ± 1.117 

2.297 ± 0.170 

14900 ± 1100 

06-04T03:35:31 

DUAL 

31.46 

12.861 ± 2.172 

1.226 ± 0.207 

13000 ± 2200 

06-04T04:16:38 

DUAL 

11.09 

6.119 ± 1.217 

1.654 ± 0.329 

16300 ± 3200 

06-04T04:16:38 

DUAL 

31.46 

11.857 ± 3.149 

1.130 ± 0.300 

11900 ± 3200 

FOF177 






06-06T04:15:14 

SWC 

19.71 

18.338 ± 2.343 

2.789 ± 0.356 

17900 ± 2300 

FOF178 






06-llT03:53:40 

SWC 

19.71 

16.962 ± 2.246 

2.580 ± 0.342 

16100 ± 2100 

FOF179 






06-13T04:04:51 

SWC 

19.71 

16.991 ± 1.777 

2.584 ± 0.270 

16000 ± 1700 


a Measured in a circular aperture with a radius of 9.984" centroided on the photocenter 
of the comet nucleus. 
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Table 3. Best-fit Thermal Model Parameters and Derived Grain Mineralogy of Comet 

C/2012 K1 (Pan-STARRS)* 


Dust component 

N p a xlO 16 

Sub-//m 
mass fraction 

Amorphous pyroxene 

8526.433™;“' 

0 8in +0 - 043 

u.oiu_ 0 060 

Amorphous olivine 

1 998 qoc;+213.563 
izzo.ozu_ 640 688 

D 04.e;+0.009 
o.uuo_ 0 026 

Amorphous carbon 

15246.5601™ 

n p;c,c;+0.009 
u.ooo_ 0 017 

Crystalline olivine 

2476.880™“' 

0 090 + °- 078 
u.uyu_ 0 031 

Crystalline pyroxene 

0.0001“ 

o 

o 

o 

o 

1 + 

o o 

b b 

O 

O CO 

Other model parameters 

xl 

0.98 


Degrees of freedom 

156 


Total submicron grain mass b 

(7.663«:i‘°) x 10 s kg 


Silicate/carbon ratio 

0-80l°;i 


fcryst 

0.202lg;SS 



*Notes. Uncertainties represent the 95% confidence level. 

a Number of grains at the peak of the grain size distribution. 

b The total mass of the sub-/im sized grains contained within the 
spectral extraction aperture. 
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Table 4. Thermal Modeling Dust Characteristics of Select Comets 


Comet Class 

1/ Cl 0 rig 

(i(r 6 Air 1 ) 

f a 

Jcryst 

Range 

(%) 

SAC b 

Ratio 

N 

&peak 

(/mr) 

Refs* 

NIC/OC c 

C/2012 K1 (Pan-STARRS) 

42.9 

10-50 

0.6-1.1 

3.4 

0.6 

This work 

C/2007 N3 (Lulin) 

32.2 

34-51 

0.42-0.54 

4.2 

0.9 

[1] 

C/2001 Q4 (NEAT) 

61.2 

71 

1.7-5.7 

3.7 

0.3 

P.3] 

C/2002 VI (NEAT) 

2279.3 

66-69 

1.38 

3.7 

0.5 

[3] 

C/1995 01 (Hale-Bopp) 

3805.0 

60-78 

8.1-13.3 

3.4-3.7 

0.2 

[4,5] 

EC/JFC C 

9P/Tempel 1 (~ lhrs post-impact, ctr) 


19-25 

3.4-4.4 

3.7 

0.2 

[7] 

17P/Holmes 


~42 c 

0.2 



[7] 

73P/SW3-B (Apert B) 


43-69 

1.09-1.59 

3.4 

0.5 

[8] 

73P/SW3-C (Apert M) 


57-69 

0.60-0.75 

3.4 

0.3 

[8] 


a Computed from thermal model sub-/im silicate mass fractions using Eqn. tITj) described in 1 13.31 
where the range includes the model uncertainties when known. 

b Dehned as the silicate-to-amorphous carbon ratio derived from thermal modelling of the SEDs 
where the range includes the model uncertainties when known. 


c Comet dynamical class divisions are NIC/OC = Nearly Isotropic/Oort Cloud; EC/JFC = 
Ecliptic/Jupiter-Family 


cl Estimated from Reach et al . (120101 ) who provide abundances weighted by grain surface area. 


References. ( TnWoodward et al.l ( 20111) ; 12) I Wooden et al.l (2004 ); (' 3)lQotsubo et al.l (2007 ): 


(41 


Harker et al 


( 2Qlol ): 18)1 Ha rker e t al. ( 2011 ) 


(120021); (5) lHarker et al.l (l2004al ): (6) lHarker et al.l (120071 ): (7) iReach et ah 
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Table 5. Afp Values for Comet C/2012 K1 (PanSTARRS) on 2014 June 04.24 UT 


Aperture a 
(arcsec) 

P 

(km) 

R 

(mag) 

Afp b 
(cm) 

11.24 

6705 

12.157i0.008 

5731 ± 42 

19.97 

11912 

11.593±0.009 

5424 ± 43 

25.14 

14999 

11.344±0.009 

5417 ± 43 

33.53 

19999 

11.065±0.009 

5253 ± 41 

41.91 

24999 

10.835±0.009 

5194 ± 41 

50.29 

29999 

10.674±0.009 

5020 ± 39 


a Effective circular aperture diameter. 

b Afp values corrected to zero phase (see H3.5|h 






